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I. Introduction 


With the availability of extensive data on two-body weak hadronic decays 
of heavy flavor mesons [1], it has now become possible to test the validity of 
the factorization model which has been considered to be supported by the D 
meson phenomenology [2-4], In the recent works [5-8], it has been shown that 
the factorization model fails to account for the observed data on charmed and 
bottom meson hadrons. Large N c limit, in which nonfactorizable contributions 
are usually ignored, does not work when extended to B meson decays as these 
clearly demand [6] a positive value of the QCD parameter 02 . Even in the D 
meson decays, universal choice of the parameters d\ and <22 does not explain 
many of the hadronic decays of D and D s mesons. For instance, the measured 
branching ratios of 77 and rf emitting Cabibbo-angle-favored decays of charmed 
mesons are considerably larger than those predicted in the spectator quark 
picture. Annihilation terms, if used to bridge the discrepancy between theory 
and experiment, require large form factors particularly for D° —> 77 / 7 / + K° 
and D° —>77 + A '* 0 decays [7]. Further, factorization also fails to relate 
Dt v/v' + tt + /p + decays with scmileptonic decays Df —> 77 / 77 ' + e + + v 
[7, 8] in a consistent manner. 

Recently, there has been a growing interest in studying the nonfactorizable 
terms for weak hadronic decays of the heavy flavor hadrons. Many attempts 
have been made to estimate the amount of nonfactrizablc effects needed to 
reproduce the experimental results for charmed and bottom sector [9-11] when 
real value N c — 3 is used. In an earlier work [12], using isospin symmetry 
one of us (RCV) has searched for a systematics in these estimates for various 
decays of D + and D° mesons. It has been shown that the nonfactorizable 
isospin 1/2 and 3/2 reduced amplitudes may bear a universal ratio for D —> 
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Kir / Kp/K*n / Ka,\/K*p decay modes. The formalism, when generalized to 
SU (3) for studying D —>• PP decays, where P denotes a pseudoscalar meson, 
has resulted in a consistent fit with experiment [13]. 

In this work, we have extended the SU (3)-flavor analysis of nonfactorizable 
contributions to D —» PV decays, where V denotes a vector meson. We analyze 
the Cabibbo- favored decays of D°, D + and D F mesons taking into account the 
final state interactions (FSI). In section II, we develop the formalism. Results 
and discussion are given in the last section. 

II. Formalism 

We start with the effective weak Hamiltonian 

H w = G F [ci(ud)(sc) ± c 2 (sd)(uc)], (1) 

where Gp = V u dV* s and q\q 2 (= — 75 )^ 2 ) represents color singlet V - 

A current and the QCD coefficients at the charm mass scale are 

ci = 1.26 ±0.04, c 2 =-0.51 ±0.05. (2) 

Due to the Fierz transformation of the product of two Dirac currents in (1) in 
W'Color space, the Hamiltonian takes the following form [9]: 

H% F = G F [ai(ud)(sc) + c 2 H*\, 

H^ 3 = G F [a 2 {sd)(uc)+ Cl H *], (3) 

for color favored (CF) and color suppressed (CS) decay amplitudes respectively. 
Here, 

01,2 = ci , 2 ± ( 4 ) 

Hi = A(sw)(a*c), 

Z a= 1 
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( 5 ) 


Hi = \ 't(sX‘d)(u\‘c), 

Z a= 1 

where qi\ a q 2 (= < 717^(1 — 75 )A a g 2 ) represents color octet current. 

Matrix elements of the first terms in (3) can be calculated using the factor¬ 
ization scheme. These are given column (ii) of Table 1. So long as one restricts 
to the color-singlet intermediate states, second terms in (3) are ignored and 
one usually treats oq and a 2 as input parameters in place of using N c = 3 in 
reality, ft is generally believed [2-4] that the D —> Kn decays favor N c —» oo 
limit, i.e., 

cq « 1.26, a- 2 ~ —0.51. ( 6 ) 

However, it has been shown that this does not explain all the decay modes of 
charm mesons [5,7]. For instance, the observed D° —» r)K 0 and D° —» r/K° 
decay widths are larger than those predicted in the spectator quark model. 
Also in D — > PV mode, measured branching ratios for D° —> K°uj / r)K *°, 
D + —> K*°n + , and Df —> K°K* + /qp + /rf + p + , are considerably higher 
than those predicted by the spectator quark diagrams. In addition to the 
spectator quark diagram, factorizablc W-exchange or W-annihilation diagrams 
may contribute to the weak nonlcptonic decays. However, such contributions 
are normally expected to be suppressed [2] in the meson decays. For D + 
meson decays, these do not appear in the Cabibbo favored decay process. For 
D° meson decays, these are further color-suppressed as these involve lower 
QCD coefficient o 2 . Therefore, we have ignored them in the present analysis. 

We now investigate nonfactorizable contributions to these decays. Matrix 
elements of and FT.® between charm mesons and two-body uncharmed final 
states are difficult to calculate theoretically [9,10], as these involve nonper- 
turbative effects arising due to soft-gluon exchange. We employ SU(3)-flavor- 
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symmetry [14] to handle these matrix elements. In the SU(3) limit, the two 
Hamiltonians and H ’® behave like Hf 3 component of 6* and 15 represen¬ 
tations of the SU(3). Since Hf i; and Hf,, transform into each other under the 
interchange of u and s quarks, which forms V-spin subgroup of the flavor- 
SU(3), the reduced matrix elements satisfy 

<PV\\Hl\\D> = <PV\\Hl\\D>. (7) 


The matrix elements < PV\Hf u \D >, appearing in the nonfactorizable effects, 
are considered as weak spurion + D meson —> P + V scattering process, 
whose general structure can be written as 


< PV\Ht\D > = h(P‘ V™P C ) + a 2 (P™ViP<) 


+«s 1 

UMahcn+ 

+Mftp") + h(P^VfP")]Hfa 
+\^P c XP b ) + 62 (PXP^Hf^ 
+\d l (p;y;;:p i ) + 


( 8 ) 


where P a denotes triplet of D-mesons P a = ( D °, D + , £H~) and P b “, 14“ denote 
3 (8) 3 matrices of uncharmed pseudoscalar meson and vector meson nonets 
respectively. For pseudoscalar mesons, 



(Pi 

7T + 

K+\ 


pa _ 

r b ~ 

71 

p2 

K° 

(8) 


\K~ 

K° 

Pi) 



with 

Pi = —j= (7r° T T) sin 6 -\- fj' cos 0 ), 
V2 
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P 2 2 = —(—7T° + rj sin 9 + r\ cos 6 ), 
v2 

P 3 3 = — rj cos 6 + 77 ' sin 9, (9) 

where 9 governs the rj — 7/ mixing, and is related to the physical mixing [1] as, 

0 = 9 idea i — 4>ph.y ( 10 ) 


For vector meson nonet, 



( v? 

P + 

K *+ \ 


V b a = 

p 

vi 

K*° 

(ii) 


\K*- 

K*° 

vi ) 



with 

Vi = —=(p° + u;cos9' + (j)sm9')), 

V2 

V? = —p(— p° + u>cos9' + 4>sm9')) i 
v2 

Vij 3 = uj sin 9' — cj) cos9'. 

For ideal to — (j> mixing 

9' = 0. (12) 

In addition to the nonfactorizable effects considered so far, there may also 
arise nonfactorizable effects involving product of color-singlet currents. How¬ 
ever, these may be relatively suppressed [10]. Even if these are included here, 
it has been shown [13] that their contributions can be absorbed in the unknown 
reduced amplitudes appearing in (7) due to the similar structure in the SU(3) 
framework. 

There exists a straight correspondence between the terms appearing in 
(7) and various quark level processes. The terms involving the coefficients 
ai, a 2 , bi and b 2 represent annihilation diagrams. Notice that, due to the 
involvement of gluons, these are no longer suppressed. The terms having 
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coefficients 03 , 63 and 64 represent spectator-quark like diagrams where the 
uncharm quark in the parent D-meson flows into one of the final state mesons. 
The last terms having coefficients e\s and d's are hair-pin diagrams, where 
qq pair generated in the process hadronizes to one of the final state mesons. 
Choosing Hf 3 component from (7), we obtain nonfactorizable contributions to 
various D —> PV decays. These are given in column (ii) of Table II. 

III. Results and Discussion 

The decay rate formula for D —> PV decays is given by 

T(D —> PV) = \G f \ 2 -^—\A{D^PV)\\ (13) 

oTTTTly 

where p c is the three-momentum of final state particles in the rest frame of 
D meson and my is the mass of vector meson emitted. Now we proceed to 
determine nonfactorizable effects to various decays. First, we determine the 
factorizable contributions to various decays using N c = 3 which fixes, 

q 1 = 1.09, Qj'2 — — 0.09, (14) 

ignoring the errors in the QCD coefficients Ci and C 2 . For the form factors at 
zero momentum transfer, we use 

F? k { 0) = 0.76, F° n ( 0) = 0.83, 

A£ K \ 0) = 0.75, (15) 

as guided by the semileptonic decays of D— mesons [1, 15-17], and 

F, Dt, ( 0) = 0.681, F, Dr/ (0) = 0.655, 

F? bT, ( 0) = 0.723, F? W ( 0) = 0.704, 

F° sK ( 0) = 0.760; (16) 
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Aq P { 0) = 0.669, Ag u ( 0) = 0.669, 

A^( 0) = 0.669, A° s “( 0) = 0.700, 

Ao s *( 0) = 0.700, A° aK \ 0) = 0.634, (17) 

from the BSW model [2] for want of experimental information. These form 
factors are then extraploated in q 2 using a monopole form with appropriate 
pole masses. Following values are used for the meson decay constants [1,2] (in 
GeV) 

V = 0.132, f K = 0.161, 

f P = 0.212, f K . = 0.221. (18) 

Numerical values of the factorizable amplitudes are given in col (iii) of Table 
I. 

Notice that D —> Kp and D —> K*n decays involve elastic FSI whereas the 
remaining decays are not affected by them. As a result, the isospin amplitudes 
1/2 and 3/2 appearing in these decays may develop different phases. We 
illustrate the procedure for D —> Kp amplitudes: 


b 

o 

l 

- K~p + ) 

= -^= [A 3/2 e^/ 2 + v / 2A 1/2 e i<5l/2 ], 


A(D° - 

K°p°) 

= -L[v^A 3/2 e^/ 2 -A 1/2 e^/ 2 ], 



A(D + - 

*K°p+) = VS A 3/2 e iS ^. 

(19) 


Following phase independent relations: 

\A{D° -> K ~ p +)\ 2 + \A{D° -> iF°p°)| 2 = |A 1/2 | 2 + |A 3/2 | 2 , 

\A(D + —>■ K ° p + )\ 2 = 3|A 3/2 | 2 , (20) 

allow us to work without the phases. Writing the total decay amplitude as a 
sum of the factorizable and nonfactorizablc parts 

AID — PV) = A f (D — PV) + A n, (D — PV), 


( 21 ) 



we obtain 


A™j 2 (D Kp) = ~^{V2A nf (D°^K-p + )-A nf (D 0 ^K°p 0 )}, (22) 


V3 


A”f 2 (D -> Kp) = —{A nf (D° K~p + ) + V2A nf (D° -*■ iP 0 p°)} 


v/3 


V3 


A n/ (D + ->■ ji'V). 


(23) 


The last relation (23) leads to the following constraint: 
AJ/ 2 (B -* tfp) c? + 2ci 


- tfp) 


\/2(c| - cl) 


= - 1.123 ±0.112. 


(24) 


Experimental value B(D + —> K°p + ) = 6.6% then predicts sum of the branch¬ 
ing ratios of H 0 —» Kp modes: 


5(£> u -> iP"p + ) + 5(L>° -> K°p°) = 10.42%, (11.50±1.31% Expt.) (25) 


in good agreement with experiment. Present data [1] on these modes is consis¬ 
tent with a choice of zero phase difference between isospin 1/2 and 3/2 channels 
and fixes the ratio of nonfactorizable amplitudes, 


A r ;jjD Kp) 

-VL-—- = - 1.481 ± 0.582 (Expt 

A" 3 j 2 (D - Kp) 


(26) 


consistent with theoretical value given in Eq. (24). For branching ratio of the 
individual modes, we obtain 


£(£>° - 

-> K~ p + ) = 9.37% 

(10.4 ± 1.3% Expt.), 

(27) 

£(D° - 

^/?V 0 ) = 1-05%, 

(1.10 ±0.18% Expt.). 

(28) 


Performing a similar analysis for the D —> K*n decay amplitudes and using 
the experiemental value B(D + —> K*°7 r + ) = 2.2%, we calculate: 

B(D° —> K*~tt + ) ± B(D°->K*°7 T°) = 7.44%, (7.9 ±0.7%, Expt.) (29) 
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in nice agreement with experiment. For these modes also, the isospin reduced 
amplitudes bear the same ratio as given in Eq.(24), 


A r/ 2 (r» -> »f) _ J + 2cl 


V2(c% - cf) 


= - 1.123 ±0.112, 


(30) 


which compares well with experimental value —1.171 ± 0.158 when negative 
sign is chosen for A$/ 2 - 

Calculation of branching ratio of the remaining D —> PV decays needs 
numerical values of the reduced amplitudes. Apparently these decays seem 
to involve several unknown parameters. However, the parameters d's and e’s 
appear only when isosinglet meson is emitted. Also that D meson decays 
involve combinations (cq — b \) and (02 — 62 ) which in fact are expressible as 


ai — 61 — [—(ci ± 02)03 ± C263 + C\b^[/(c\ — C2), 

0-2 — 62 = [±(ci ± €2)03 ± C1&3 ± C264]/(ci — C2), ( 31 ) 


where 

a 3 = - 0.042 GeV 2 , (32) 

63 ± b A = - 0.251 GeV 2 , (33) 

are given by D + modes. The relations (31) follow from the constraints given in 
Eqs. (24) and (30). With the experimental values B(D° —■» K°(j)) = 0.83%, 
and B(D° —> K°co ) = 2.0%,, and taking negative and positive signs for their 
experimental amplitudes, we find (in GeV 2 ), 


63 = 0.042, 

(34) 

64 = -0.293, 

(35) 

ei + di = — 0.047. 

(36) 
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The relations given in (31) then yield: 


di — b\ — — 0.202, 

(37) 

a 2 — b-2 = 0.108. 

(38) 

Using the measured branching ratios B(Df —> K + K*°) 

K°K*+) = 4.2%, we find (in GeV 2 ), 

— 3.3% and B(D 7 —+ 

ai + bi = 0.190, 

(39) 

+ ^2 = 0.123, 

(40) 


by taking negative and positive signs of their experimental amplitudes respec¬ 
tively. These parameters then predict 

B{Df -> vr°p + ) = B{D+ 7r+p°) = 0.10% (41) 

which is well below the experimental upper limit (< 0.28%) for B(Df —> 

7r + p°). Further, experimental value B(D —> 7r + 0) = 3.5% yields, 

ei - di = 0.070 Gel/ 2 , (42) 

for the positive choice of its experimental amplitude. This value in turn leads 
to: 

B(D^n + u) = 0.61% (43) 

obeying the exprimental upper limit (< 1.7%). Now, we are left with rj 

and r)' emitting decays which invlove mixing angle (j)phy We have chosen 
to present results for all the three mixing angles —10°, — 19°, and —23° 

given in the Particle Data Group [1] so as to make the trend with mixing 
angle evident. Measured branching ratios B(D° —> rjK* 0 ) = 1.9% and 

B(D+ —> r]'p + ) = 12.0% fix the parameters: 

e 2 + d- 2 = 0.668 GeV 2 for <pph. y = — 10° 
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(44) 


= 0.373 GeV 2 for d Phy = - 19°, 

= 0.313 GeV 2 for (j> Phy = - 23°, 

and 

e 2 — d 2 = 0.457 GeV 2 for (p P h y = — 10° 

= 0.428 GeV 2 for (f> P h y = — 19°, 

= 0.412 GeV 2 for <j> Phy = - 23°, (45) 

for negative and positive signs of respective experimental amplitudes. Calcu¬ 
lated branching ratios for rj and rf emitting decays are listed in columns (ii) to 
(iv) of Table III for different values of the mixing angles. For the sake of com¬ 
parison with factorizable part, nonfactorizablc contributions to various decays 
are given in column (iii) of Table II. 
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Table I 


Spectator-quark decay amplitudes 


Amplitude (A f ) 


Process 


a i — 1.09, cl 2 — —0.09 


D° - 
D° 
D° 
D° 
D° - 
D° - 
D° 
D° - 


> K~p+ 
■> K°p° 

K°u 
-> K°(f) 
n+K*- 

> 

■» rjK*° 

> rj'K * 0 


D+ -»• K°p+ 
D+ _> tt+J?* 0 


^ s + - 

£ s + - 

^ s + 

Dt~ 

Df- 


-> 7 T ' p 
-> 7 T°p + 

-> 7r + cu 
-> 7T + 0 
-» r/p + 

-> r/p + 
K + K*° 
K°K*+ 


2a 1 f p m p Ff ,K (m 2 p ) 
\/2a 2 f K m p AQ p (m 2 K ) 
\f2a 2 f kcosQ' rn^A^ (rn 2 K ) 
-V2a 2 sin9' f K m <j> A^( m 2 K 
2a 1 f 7r m K *A^ K * (ml) 
V2a 2 f K *m K * F?*(m 2 K .) 
\/2 a 2 f K *sin9m k* F[* v (m 2 K * 
V2a 2 fx* cosOrn k* F[ Ji ' (m\ 


I\* 


?DK ( m 2) 


2a l f p m p F p 
+2 a 2 f K m p Ao P (m 2 K) 
2a 1 f n m K *A{ 2K * (ml) 
+2a 2 f K *m K *Ff >w (m 2 K ,) 

0 

0 

2a 1 f 7T sin0'm u Ao sU} (ml) 
2a l f 7T cos9'm <t> A® a(t> (ml) 
—2ai f p cos9m p Fi s11 (rn p 
2aif p sin9m p Fi sV (m 2 p ) 
2 a 2 f K * m K * F® s K ( m 2 K ,) 
2a 2 f K m K *AQ sK * (m 2 K ) 


0.316 

- 0.011 

- 0.011 

0 

0.179 

-0.026 

-0.015 

-0.015 


0.300 

0.143 

0 

0 

0 

0.204 

-0.209 

0.205 

-0.033 

-0.018 
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Table II 


Nonfactorizable contributions to D —> PV decays, 
(ci and c 2 are the QCD coefficients) 


Process 

Nonfactorizable contribution 


D° - 

4 K~p+ 

C2[-ai - a 3 + bi + 64] 

0.024 

D° - 

4 K°p° 

-^2 ci [a 1 -a 3 -bi + b 3 } 

- 0.103 

D° - 

-4 K°to 

ci[ ( a i a 3 “t bi T ^3 4 " 2 (ei + di)) + sinO ( 02 + b 2 + ei + d±)] 

0.174 

D° - 

-4 K°<t> 

Cl [y^-(- ai - «3 + h + 63 + 2(ei + di)) - cos 0 '{a 2 -b 2 - e 1 - di)] 

- 0.199 

D° -4 

■ n + K*~ 

C2[-«2 + a 3 + 62 + 63] 

0.058 

D° - 

> tt°K*° 

■^jCi[a 2 + a 3 — 62 + ^4] 

-0.201 

D° - 

4 r]K*° 

ci[^(-°2 + Q 3 + ^2 + ^4 + 2(e 2 + d2)) + co S d(ai — b\ — e 2 — c?2)] 

-0.210 

D° - 

4 n'K*° 

Ci [^^(~ a 2 + o 3 + b-2 + 64 + 2(e 2 + d 2 )) + sin 0 (—ai + 61 + e2 + c?2)] 

1.334 

D+ - 

-4 K°p+ 

(ci + C2) [ 2 u 3 + b 3 + 64] 

-0.122 

D+ - 

4 n+K * 0 

(ci + c 2 )[2a 3 + 6 3 + 64] 

- 0.254 

Df- 

-4 7 T + p° 

-^2 C2 [ a 3 + a 2 — bi + 6 2 ] 

0.024 

Df- 

-4 7 T° p + 

^C2[ai — a 2 + bi — 6 2 ] 

0.024 

Df 

- 4 - TT + U) 

c^Tif'I 01 + °2 + bi + b 2 - 2(ei - di)) + sin 0 '(a 3 + b 3 - e 1 + di)] 

- 0.062 

Df 

-4 n + (f> 

C2[ £ ^(ai + a 2 + bi + 6 2 - 2(ei - di)) - cos 0 '(-a 3 - b 3 + e 1 - di)] 

0.037 

Df 

-» VP + 

c 2 ( a i + 02 + &i + b 2 — 2 (e 2 — d 2 )) + cos 0 (a 3 — 64 + e2 — d2)] 

- 0.097 

Df- 

~4 V’P + 

c 2 [ c ^ (ci + ci2 + 61 + b 2 — 2(e 2 — d 2 f) + sinQ( — a 3 + 64 — e 2 + d 2 )\ 

0.424 

Df - 

■ K+K*° 

ci jai + a 3 + b\ + 64] 

- 0.186 

Df - 

■ K°K*+ 

ci [a 2 - a 3 + b 2 + b 3 } 

0.263 
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Table III 

Branching ratios (%) of r)/rj' emitting decays for different mixing angles. 


Decay 

<j> phy = -10° 

<l>phv = - 19 ° 

</) P hy = - 23 ° 

Experiment 

D° -> iiK*° 

1 . 90 “ 

1 . 90 “ 

1 . 90 “ 

1.9 ± 0.4 

D° -> r]'K*° 

0.33 

0.10 

0.08 

< 0.11 

D+ -► i W + 

10.47 

5.23 

3.50 

10.0 ± 2.2 

D+ -» rfp+ 

12 .00“ 

12 .00“ 

12 .00“ 

12.0 ± 3.0 


a input 
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